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Bovine coronaviruses (BCoVs) are widespread around the world and cause enteric or respiratory infections 
among cattle. The current study includes 13 samples from BCoVs collected in Normandy during an 11 -year period 
(from 2003 to 2014), 16 French HCoV-OC43s, and 113 BCoVs or BCoVs-like sequence data derived from partial or 
complete genome sequences available on GenBank. According to a genotyping method developed previously for 
HCoV-OC43, BCoVs and BCoVs-like are distributed on three main sub-clusters named C,, C 2 , and C 3 . Sub-cluster 
C, includes BCoVs and BCoVs-like from America and Asia. Sub-cluster C 2 includes BCoVs from Europe. Sub-cluster 
C 3 includes prototype, vaccine, or attenuated BCoV strains. The phylogenetic analyses revealed the monophyletic 
status of the BCoVs from France reported here for the first time. Moreover, BCoV exhibits a relative genetic 
stability when compared to HCoV-OC43 we previously described from the same region. The numerous recombi¬ 
nation detected between HCoV-OC43 were much less frequent for BCoV. The analysis points thus to the influence 
of different evolutive constraints in these two close groups. 

© 2016 Elsevier B.V. All rights reserved. 


1. Introduction 

Coronaviruses are the largest enveloped single-strand RNA viruses, 
ranging from 26 to 31 kilobases in size. They belong to the Coronaviridae 
family in the Nidovirales order, also including the families of Arteriviridae, 
Roniviridae, and the last described Mesoniviridae (Lauber et al„ 2012; 
Gorbalenya et al„ 2006). 

Coronaviruses are divided into four genera named Alpha-, Beta-, 
Gamma- and Deltacoronavirus based on phylogenetic distance of highly 
conserved domains (Adams and Carstens, 2012; Woo et al„ 2012). 
Betacoronavirus genus is divided into the four clades A to D. 
Coronaviruses infect a wide range of avian or mammalian species, and 
are responsible for enteric or respiratory infections (Woo et al., 2009). 
Interest in coronaviruses increased by the emergence of two human 
coronaviruses, Severe acute respiratory syndrome-related coronavirus 
(SARS-CoV) and Middle-East respiratory syndrome coronavirus 
(MERS-CoV) in 2002 and 2012, respectively. Both these coronaviruses 
belong to the Betacoronavirus genus, clades B and C, respectively, and 
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cause severe respiratory diseases in humans. SARS-CoV emerged in 
human in 2002-2003, causing an outbreak with more than 8000 cases 
with a fatality rate of nearly 10% (Peiris et al., 2004). MERS-CoV emerged 
in 2012 in the Middle-East (Zaki et al., 2012). To date, 1413 laboratory 
cases of MERS-CoV infection have been confirmed, including 502 fatal 
cases, reported by the World Health Organization (WHO) (Institut 
National de Veille Sanitaire (INVS), 2015). Both SARS-CoV and MERS- 
CoV have a zoonotic origin, with an animal reservoir in close contact 
with the human population. The zoonotic potential of coronaviruses 
implies the need for surveillance of coronaviruses associated with 
domestic animals in close contact with human population. 

Bovine coronavirus (BCoV) belongs to the Betacoronavirus genus 
clade A. This clade also contains other closely related coronaviruses 
known as Bovine-like coronavirus (BCoV-like), originating from captive 
wild ruminants such as waterbuck (Kobus ellypsiprimus), sambar deer 
(Cervus unicolor), white-tailed deer ( Odocoileus virginianus), elk ( Cervus 
elephus), giraffe ( Giraffa Camelopardalis), and sable antelope ( Hipotragus 
niger) (Alekseev et al„ 2008). These coronaviruses were antigenically, 
biologically, and genetically highly related to BCoV (Alekseev et al., 
2008; Tsunemitsu et al., 1995; Majhdi et al., 1997; Hasoksuz et al., 
2008). BCoVs-like were also detected from camelids such as alpaca 
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(Vicugna pacos ), llama ( Lama glama ), and dromedary camel ( Camelius 
dromedarius). Particularly, a BCoV-like, named Dromedary camel corona- 
virus (DcCoV), was detected from dromedary camel feces in the United 
Arab Emirates in 2013, during an epidemiological investigation aiming 
to identify the animal reservoir of MERS-CoV (Cebra et al„ 2003; Woo 
et al„ 2014). BCoVs-like are globally responsible for enteric infection. 
However, in 2002 a BCoV-like responsible for respiratory infection was 
detected from a dog in United Kingdom (Erles et al„ 2003). To date, 
only 15 complete genomes of BCoVs, including 13 from USA, one from 
Canada, and one from South Korea are available on GenBank. Moreover, 
16 BCoVs-like complete genomes, including 12 from captive ruminants, 
one from a dog, and three from dromedary camels are available on 
GenBank (consulted on September 8 th 2015). Much earlier, in 1988, a 
coronavirus closely related to BCoV was isolated from the feces of a 
diarrheic child in Germany. This coronavirus, tentatively named at first 
Human enteric coronavirus (HEC 4408), revealed itself to be closer to 
BCoVs than to BCoV-Like and to other HCoVs indeed. The phylogenetic 
position of this virus and the context of this case suggests the infection 
of the child by a cattle virus, consequently to a species barrier-crossing 
(Zhang etal., 1994). 

BCoV is closely related to the Human coronavirus OC43 (HCoV- 
OC43), that was first isolated in 1967 (McIntosh et al., 1967). More 
precisely, BCoV and HCoV-OC43 share a global nucleotide identity 
of 96%. The genetic comparison of BCoV and HCoV-OC43 reveals 
that HCoV-OC43 results from a zoonotic transmission from bovine to 
human. Using a molecular clock analysis, Vijgen et al. estimated 
the tMRCA of these viruses around 1890 (Vijgen et al„ 2005). The intra¬ 
specific variability of the HCoV-OC43 revealed the co-circulation of 
numerous recombinant variants (Kin et al., 2015). Given the similarity 
of HCoV-OC43 with BCoV, we were interested in testing the genetic 
variability of BCoV using the same methods. Thus, the objective of the 
current study was to investigate the intraspecific genetic variability of 
BCoV collected in Normandy from 2003 to 2014 and to test whether 
recombination remain a major mechanism of variability in BCoVs as 
seen for his sister clade HCoV-OC43. 


2. Materials and methods 

2.1. Samples and positive control 

Thirteen field samples positive for BCoV and one cell-culture BCoV 
supernatant were included in this study. The 13 field samples were 
provided to us by laboratories of Normandy specializing in veterinary 
diagnostics. These samples include 12 diarrheic fecal samples and one 
nasal aspirate. The cattle herd of origin of nine samples are known pre¬ 
cisely and depicted in Fig. 1. The cell-culture supernatant is derived from 
a fecal sample propagated on HRT-18 cells in our laboratory since 1990s. 
We used this BCoV as positive control. Age, type of samples and counties 
of origin of specimens, when available, are indicated in Table 1. 


2.2. Sample preparation, nucleic-acid extraction, and real-time RT-PCR 

Fecal samples were resuspended (10% suspension in sterile PBS 
Buffer) and centrifuged at 3000 rpm for 10 min. Total nucleic acids of 
200 pi of fecal suspension, nasal aspirate, and culture supernatant 
were manually extracted using High Pure RNA Isolation Kit (Roche, 
Basel, Sweden) following the manufacturer's instructions. The determi¬ 
nation of the cycle threshold (Ct) was performed using an in-house real¬ 
time RT-PCR, targeting a 128 nucleotide fragment of the BCoV M gene, 
for which the original primers and probe can be found in Table SI 
(supplementary material). RT-PCR reactions were performed using 
One Step RT-PCR System (Qiagen, Holden, Germany). Reaction mixtures 
were performed in a final volume of 25 pi, composed of 5 pi of 5 x reac¬ 
tion buffer, 1 pi of a dNTP solution (0.4 mM each), 1 pi of Enzyme mix 
(Reverse transcriptase and Taq polymerase), and 12 pi of RNA free 
water, provided in the kit. We added 1.5 pi of MgS0 4 , 0.4 pM of each 
BCoV-M550F and BCoV-M650R primers, 0.2 pM of BC 0 V-M 6 OO probe, 
and 2.5 til of RNA template. Thermal cycling involved 50 °C for 30 min, 
followed by 95 °C for 15 min and then 45 of 95 °C for 30 s and 60 °C 
for 1 min. 



Fig-1- Location of the livestock where nine of our samples were collected. The name of sample was indicated above each farm of origin as follows: Caen/year of sampling/specimen number. 
The nine corresponding farms were named FI to F9. 
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Table 1 

Features of the 14 BCoV samples used in this study. a mo„ month; d„ days. b na, not 
available. 


Specimen 

Age at time 
sampling a,b 

Type of sample 

Department of 
origin 15 

Caen/2005/01 

na 

Respiratory 

Calvados 

Caen/2005/02 

na 

Enteric 

Calvados 

Caen/unknown/03* 

na 

Culture supernatant (HRT18) 

na 

Caen/2008/04 

<1 mo. 

Enteric 

Manche 

Caen/2003/05 

na 

Enteric 

Calvados 

Caen/2010/06 

5 d. 

Enteric 

Manche 

Caen/2012/07 

<1 mo. 

Enteric 

Manche 

Caen/2013/08 

15 d. 

Enteric 

Manche 

Caen/2013/09 

<1 mo. 

Enteric 

Manche 

Caen/2013/10 

<1 mo. 

Enteric 

Manche 

Caen/2013/11 

<1 mo. 

Enteric 

Manche 

Caen/2014/12 

<1 mo. 

Enteric 

Manche 

Caen/2014/13 

<1 mo. 

Enteric 

Manche 

Caen/2004/14 

7 d. 

Enteric 

Calvados 


* Propagated on cell culture since 1990s. 


2.3. Sequencing of Che three complete genes nsp 12, S, and N and genotyping 
method 

The sequences of BCoVs were named according to the following 
nomenclature: Virus name/FRA-EPI/location of sampling/year of 
sampling/specimen number. FRA corresponds to the acronym of 
France, EPI is an abbreviation for the EPICOREM consortium. In this 
study, the term “molecular isolate” was used for sequences of BCoV 
amplified directly from samples. 

To determine the genomic sequence of nspl2, S, and, N genes, a set 
of 17 overlapping RT-PCR products, from 553 to 1260 nucleotides, were 
generated using a total of 42 primers (Lau et al., 2011; Hasoksuz et al., 
2002; Martinez et al., 2012). RT-PCR reactions were performed using 
One Step RT-PCR System (Qiagen, Holden, Germany), from 2.5 pi of 
nucleic acid in a final volume of 25 pi, according the manufacturer's 
instructions. Briefly, reaction mixtures were composed of 5 pi of 5 x 
reaction buffer, 1 pi of a dNTP solution (0.4 mM each), 1 pi of Enzyme 
mix (Reverse transcriptase and Taq polymerase), and 13.5 pi of RNA 
free water, 0.4 pM of forward and reverse primers (Table SI, supple¬ 
mentary material), and 2.5 pi of RNA template. Thermal cycling involved 
50 °C for 30 min, followed by 95 °C for 15 min and then 45 of 95 °C for 
30 s, 56 °C to 60 °C (depending of primers used), and 72 °C for 1 min. 

The RT-PCR products encompassed all three complete genes 
(Table SI, supplementary material). All the primers were tested in silico 
and in vitro with our control BCoV strain to test their efficacy and the 
primer abilities to provide us with quality sequencing products. 

The bidirectional sequencing was performed at the Laboratory for 
Urgent Response to Biological Threats, Environment and Infectious 
Risks Unit, Institut Pasteur (Paris, France), following their sequencing 
protocol. The analysis of sequencing products was performed by the 
capillary electrophoresis 3730XL (Life Technologies, Carlsbad, CA, US). 


2.4. Bioinformatic analysis 

Sequences were assembled in contigs corresponding to the entire 
nspl2, S, or N genes with CodonCode Aligner Software, version 5.0.1 
(CodonCode corporation, Centerville, MA, USA). Multiple sequence 
alignments were performed using Muscle algorithm, implemented in 
MEGA software, version 6.06 (http://www.megasoftware.net). Thirteen 
BCoV sequences from the USA, available in Genbank, were added to the 
alignment (Tables S2 and S3, supplementary material) (Chouljenko 
et al., 2001; Zhang et al., 2007; Mebus et al., 1973). Two additional 
prototype BCoV sequences from Japan and Canada, named Kakegawa 
and Quebec respectively, were also added (Gelinas et al., 2001; Yoo 
and Pei, 2001). Sixteen sequences from bovine-like coronaviruses 
were added as well, corresponding to coronaviruses collected from 


sable antelope, white-tailed deer, sambar deer, waterbuck, giraffe, 
alpaca, dog, and dromedary camel (Alekseev et al., 2008; Cebra et al., 
2003; Woo et al„ 2014; Erles et al., 2003; Hasoksuz et al., 2007). Finally, 
the sequence of HEC 4408 published by Zhang et al. and the 16 
sequences of HCoV-OC43 previously published by our team were 
added (Zhang et al., 1994; Kin et al., 2015). The accession numbers 
corresponding to these sequences are given in Tables S2 and S3 (supple¬ 
mentary material). Note that there are no complete genomes of BCoV or 
BCoV-like isolates from Europe. 

Maximum likelihood phylogenetic trees were constructed using the 
Tamura-Nei’s (TN93) substitution model, according to Akaike Informa¬ 
tion Criterion (AIC) and Bayesian Information Criterion (BIC) values 
evaluated by a model testing (Aldrich, 1997; Tamura and Nei, 1993). A 
bootstrap test including 1000 replicates was performed for each tree 
(Felsenstein, 1985). All of these programs are implemented in MEGA6 
software (Tamura et al., 2013). The topologies of trees were also 
double-checked using a neighbor-joining method (substitution model 
TN93), with a bootstrap test including 1000 replicates (Felsenstein, 
1985; Saitou and Nei, 1987). The topologies of these trees constitute 
the basis of our genotyping methodology. The genotype of each BCoV 
or BCoV-like is defined based on the sub-cluster in which the three 
nspl2, S and N sequences of a same BCoV (or BCoV-like) is located on 
each tree. Finally, the genotype name is constructed by the juxtaposition 
of the sub-cluster names, C,, C 2 or C 3 , on nspl2, S, and N tree 
respectively. 

An additional maximum likelihood tree based on 143 complete 
sequences of S gene of BCoV, BCoV-like, and HCoV-OC43 from USA, 
Brazil, South Korea, China, Italy, Swedish, Denmark and France was con¬ 
structed. To this end, 81S gene sequences from BCoVs and BCoVs-like 
were added to the previous S gene alignment (Martinez et al., 2012; 
Beaudeau et al., 2010; Decaro et al., 2009a; Jeong et al., 2005; Zhang 
et al., 1991; Park et al., 2007; Chung et al., 2011; Park et al., 2006). The 
accession numbers corresponding to these sequences are provided in 
Tables S2 and S3 (supplementary material). Phylogenetic trees were 
constructed using the maximum likelihood method with the substitu¬ 
tion model TN93, implemented in MEGA6 (Aldrich, 1997; Tamura and 
Nei, 1993; Tamura et al., 2013). The topology of the tree was checked 
using a neighbor-joining (substitution model TN93) method, with a 
bootstrap test including 1000 replicates (Felsenstein, 1985). 

3. Results 

The phylogenetic analysis presented in Fig. 2 shows the nspl2, S, and 
N trees obtained by a maximum likelihood method. Two clusters, one 
grouping BCoVs and BCoVs-like and the other one containing HCoV- 
OC43s, are observed in the nspl2 and N trees. DcCoVs differentiates as 
a potential third cluster on the S tree. The BCoV cluster of the nspl2, S, 
and N maximum likelihood trees is divided into three sub-clusters. 
These sub-clusters are globally comprised of by the same molecular iso¬ 
lates for nspl2, S, and N trees. Nspl2, S, and N produced comparable 
tree topologies and therefore a globally congruent phylogenetic signal 
that does not suggest BCoV recombination events, contrary to HCoV- 
OC43 during the same period and in the same geographic area (Kin 
et al., 2015). The positions of BCoV and BCoV-like sequences on each 
tree are summarized in Table S2 (supplementary material). According 
to nspl2, S, and N trees, one sub-cluster named Q contains BCoVs and 
BCoVs-like from the USA, constituting the genotype QCiCi. A second 
sub-cluster, named C 2 , groups on each trees all the BCoV sequences 
generated in this study defining the genotype C 2 C 2 C 2 . Finally, the last 
sub-cluster, named C 3 , contains the prototype strains Kakegawa, 
Mebus, Quebec, and CrCoV on all trees, constituting the genotypes 
C 3 C 3 C 3 . Two exceptions are observed. First, the three DcCoVs are includ¬ 
ed on sub-cluster C 2 on nspl2 and N tree, but are grouped in a potential 
outlier position from BCoV and HCoV-OC43 clusters, constituting a 
genotype named C 2 C 0 utuerC 2 . The second exception is the HEC 
4408 that is located on sub-cluster C 3 on nspl2 and N trees, and in 
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Fig. 2. Phylogenetic analysis of the complete nspl2, S, and N genes of 29 BCoVs, 16 HCoV-OC43s, 17 bovine-like CoVs. The phylogenetic trees were constructed by the maximum likelihood 
method. Bootstrap values were calculated from 1000 replicates. Bootstrap values over 70% are shown. The evolutionary distances were computed using the Tamura Nei model and units 
are the number of base substitutions per site. 


sub-clusters C 2 on S tree, constituting the genotype C 3 C 2 C 3 . This analysis 
reveals a potential geographical distribution with a segregation of BCoV 
and BCoV-like sequences from US on one hand, and BCoV from France in 
another hand. An additional S tree, depicted in Fig. S4 (supplementary 
material), includes additional BCoV and BCoV-like sequences from the 
USA, Brazil, South Korea, China, Sweden, Denmark, and Italy. The topol¬ 
ogy of this additional S tree is concordant with the S tree presented in 
Fig. 2. The same tree clusters (BCoV and BCoV-like, HCoV-OC43, and 
DcCoV) were observed. In the same way, the three sub-clusters C lt C 2 , 
and C 3 were observed inside to the BCoV cluster. The additional 
sequences include in this most detailed S tree allowed us to deeply 
characterize the geographic distribution of BCoV and BCoV-like. Indeed, 
in addition to USA sequences, the sub-cluser C 3 contains BCoVs and 
BCoVs-like from both America and Asia continent, more precisely 
from China, South Korea, and Brazil. The sub-cluster C 2 contains in addi¬ 
tion to the 14 French BCoVs, 35 additional BCoVs from other European 
countries such as Sweden, Denmark, and Italy ( cf. Tables S2 and S3 in 
supplementary material). 

4. Discussion 

The world population of cattle is estimated at 1.43 billion (Robinson 
et al„ 2014). In France, the cattle population extends to around 19 
million among which 1.9 million constitute the livestock of Lower 


Normandy, where the present study emerges from. BCoV infections 
are not listed as “monitored transmittable veterinary infections". There¬ 
fore, no public data measuring the impact of this infection on livestock 
cattle, national or international, are available despite BCoV being an 
ubiquitous virus. Partial genomic sequences available in GenBank 
come from 17 countries, unequally distributed on five continents. 
Fifteen complete BCoV genomes from North America and South Korea 
(14 and 1, respectively) are currently available (according to GenBank, 
consulted on September 8th, 2015) but no BCoV sequence from 
France was available before the present study. Among the few available 
epidemiological studies, Beaudeau et al. showed that in Sweden, the 
seroprevalence of BCoV is ranging from 43 to 65% in cattle less thanl2 
months, and observed a winter circulation in the temperate zones of 
the Northern hemisphere (Beaudeau et al., 2010). This epidemiology 
shares common characteristics with that of the HCoV-OC43 infection 
in humans, namely a primary infection in young patients in winter. 
BCoV and HCoV-OC43 belong to the clade A of Betacoronavirus genus 
sharing a global similarity of 96% on a nucleotide level. The divergence 
between HCoV-OC43 and BCoV is estimated at the end of the 19th cen¬ 
tury, thus contemporary with a species barrier crossing from cattle to 
humans. Six different genotypes have been identified in France among 
the 15 HCoV-OC43 circulating strains (including numerous recombi¬ 
nant variants) in respiratory samples taken over a 13-year period (Kin 
et al., 2015). This evolutive pattern (deeply affected by recombination) 
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has been commonly described for other coronaviruses ( Infection 
bronchitis virus, Feline coronavirus, Canine coronavirus ) and proves, in 
these contexts, to be a factor that favors the potential emergence of 
new variants (Cavanagh et al., 1992; Herrewegh et al„ 1998; Decaro 
et al., 2009b). There is no published data for the potential localization 
of recombination hot spot(s) on the BCoV genome. Given the genetic re¬ 
latedness between HCoV-OC43 and BCoV, we propose the hypothesis 
that these potential hot spots could be located in the same regions for 
these two viruses. Extensive sequencing of complete BCoV genomes 
would ascertain the existence of these hot spots, but was not performed 
in this study. 

The addition of BCoV sequences in France into the phylogenetic 
analysis does not modify global phylogeny since HCoV-OC43 and 
BCoV remain monophyletic in our analysis. This holds true even when 
the prototype strains (propagated in cellular culture) are included in 
our analysis (data not shown). In this light, BCoVs from France are 
grouped in a non-ambiguous way with the other BCoVs but remain 
themselves monophyletic. Base on complete or partial genome phylo¬ 
genetic trees, BCoVs-like form a sister group to BCoVs and include 
BCoVs from various animal origins (Alekseev et al., 2008; Woo et al., 
2014). 

The topological analysis of the three phylogenetic trees obtained, 
nspl2, S, and N, shows the existence of two clusters, one including 
HCoV-OC43 and the other including both BCoV and BCoV-like, respec¬ 
tively. The BCoV cluster is divided into three sub-clusters: C,, C 2 , and 
C 3 . Results obtained using the chosen methodology do not show the 
existence of circulating BCoV recombinant variants. Indeed, all of the 
BCoVs in our study are of non-recombinant genotype C 2 C 2 C 2 (see 
Table 2). All of the sequences of the 24 BCoVs and BCoVs-like coming 
from the US and included in the phylogenetic analysis are of a non¬ 
recombinant CtC^! genotype. In the same way, the sequences of 
three BCoV prototypes and the CrCoV included in the study are of a 
non-recombinant C 3 C 3 C 3 genotype. Within this analysis of putative 
recombinant detection, only the coronavirus strain HEC 4408 shows a 
recombinant genotype of type C 3 C 2 C 3 . It is worthwhile to highlight 
that the HEC 4408 finds his root among the French isolates of BCoV 
(Fig. 2) and thus represents an independent sub-group within the 
BCoV cluster. The HEC 4408, serologically distinct from HCoV-OC43, 
illustrates the capacity of BCoV for spill-over independently of its role 
in the emergence of HCoV-OC43. 

In the end, these results suggest that HCoV-OC43 and BCoV do not 
seem to evolve along the same dynamics, as it points to a generation 
of recombinant variants that seems very active in HCoV-OC43 as 
compared to BCoV. Given that these two viruses are extremely close 
with respect to genetic and epidemiological characteristics, this 
difference could be linked to different environmental constraints 
( sensu-lato ). In this light, one may conceive of a lesser genetic diversity 
of the animal host (animal livestock versus human) resulting from 
domestication and breeding practices by human, and a lesser geographic 
mobility of animals, thereby weakening the BCoV diversifying selection 
process and decreasing potential co-infection and thus recombination 
events. These hypotheses must however be tested on a greater viral 
genetic dataset. 

The results involving DcCoV sequences available in GenBank show a 
potential genotype C 2 « Coutlier » C 2 The potential outlier topology 


Table 2 

Distribution of sequences of BCoV and BCoV-like among the sub-clusters on nspl2, S, and 
N tree. 



nspl2 

s 

N 

Conclusion 

14 BCoVs from France (this study) 

Three Prototype strains* and one CrCoV 

C 2 

C 2 

C 2 

c 2 c 2 c 2 

c 3 

C 3 

C 3 

C 3 C 3 C 3 

12 BCoV and 12 BCoV-like from USA 

C, 

C, 

C, 

C,c,c, 

HEC 4408 

C 3 

C 2 

C 3 

c 3 c 2 c 3 

Three DcCoV 

c 2 

Coutlier 

C 2 

C 2 C 0 UTLIERC 2 


* BCoV Mebus, Kakegawa and Quebec. 


describes the fact that, in the S gene tree, the DcCoV sequences form a 
potential third cluster that is distinct from that of BCoV and HCoV-OC43. 

The complementary phylogenetic analysis involving a great number 
of S gene sequences presented on Fig. S4 (supplementary material) 
highlights three clusters that group each of the HCoV-OC43s, BCoVs 
and BCoVs-like, and DcCoVs, respectively. In the cluster of BCoV and 
BCoV-like strains, a geographic distribution may be noted, having a so- 
called European sub-cluster and a so-called America-Asia sub-cluster. 
This distribution suggests that the divergence of the two sub-clusters 
of BCoVs is a relatively recent phenomenon with regards to the history 
of livestock cattle. Indeed, the Bos taurus taurus livestock cattle of the 
Old World (Europe, Asia, and Africa) are issued from the domestication 
facilities of the aurochs, a wild species living exclusively in the Old 
World. This domestication event occurred around 10,000 years ago. 
Livestock cattle arrived on the North and South American continent 
after the discovery of the New World in the Xlth century (Bollongino 
et al., 2012; Loftus et al., 1994). The distinction of genetic clusters corre¬ 
sponding to the BCoV of different geographic origins could be attributed 
to the diffusion of strains favoring the exchange of animals incurred by 
international commerce. In Europe, France is the largest producer and 
exporter of cattle. Moreover, France exports live cattle for reproduction 
to Italy, Spain, Greece, and Germany, and also to Algeria, Tunisia, and 
Libya. The import-export of fresh and frozen beef between the US and 
the European Union has been minimal, especially after 1998 at which 
point the EU appealed to the World Trade Organization (WTO) to pro¬ 
hibit the import of beef for purposes of food safety. The final decision 
of the WTO was favorable to the European Union, but was not handed 
down until 2007. However, between 1998 and 2007, commerce 
between US and EU was greatly weakened. 

The United States exports beef to Asia, especially to Japan and to 
South Korea. We recall that the BCoVs, as well as the BCoVs-like collected 
from ruminants in the USA and in South Korea are very similar. It is also 
worthwhile to note that the BCoV is not included in the WTO list of 
illness being monitored in international exchange. 

In conclusion, the present epidemiological study has made it possi¬ 
ble for the first time to obtain partial sequences of BCoV originating 
from France. The results obtained, despite their limits in methodology 
and in the number of sequences, suggest the existence of an evolutive 
pattern by recombination that is different between two very similar 
human and bovine coronaviruses. This study also draws light to the 
scarce availability of BCoV sequences and genetic analyses despite the 
fact that this pathogenic agent was identified over 50 years ago and 
that it seems to have a significant sanitary impact on livestock cattle. 

Supplementary data to this article can be found online at http://dx. 
doi.org/10.1016/j.meegid.2016.03.006. 
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